MEAN-FIELD- AND CLASSICAL LIMIT OF MANY-BODY 
SCHRODINGER DYNAMICS FOR BOSONS 



Jiirg Frohlich^, Sandro Graffi^ , Simon Schwarz^ 
Abstract 

We present a new proof of the convergence of the A^—particle Schrodinger dynamics 
for bosons towards the dynamics generated by the Hartree equation in the mean-field 
Hmit. For a restricted class of two-body interactions, we obtain convergence estimates 
uniform in h, up to an exponentially small remainder. For H = 0, the classical dynamics 
in the mean-field limit is given by the Vlasov equation. 



1 Introduction and statement of results 

Consider the Schrodinger operator 



Hn = H% + Wn (1.1) 

< = -Ey^- WN = -Y.w(^i-^j) (1-2) 

i=l i<j 



where w is a two-body potential independent of N. The operator H]\f acts on TC^^^ := 
L|(M^^), the totally symmetric part of L^(R'^^), which is the Hilbert space of pure state 
vectors for a system of N nonrelativistic bosons. We propose to study the dynamics de- 
scribed by the A^— body Schrodinger equation 

ihdt^N{t)=HN^Nit), (1.3) 

for an inital condition "if]\f(t = 0) = ^'at^o S L|.(M^^). Under assumptions specified below, 



Hn, defined on the symmetrized Sobolev space H^iM.^^), is a self-adjoint operator. Hence 
the unitary group UN(t) = e"*-^^*/^, t £ R, exists. Let p< N, and let a^P) be a bounded 
operator on L|(M'^p). It defines an operator A^^^ acting on Ti^^^ in the following way: 

*(xi,...,x^)g4(m3^), (1.4) 

where Ps is the projection onto the symmetric subspace L|(IR^^) of L^(IR^^). The operator 

may be viewed as an operator acting on p particles; the numerator on the right side 
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of (1.4) is a combinatorial factor motivated by "second quantization"; the denominator is 
the correct scahng factor to take the N ^ oo hmit. 

We are interested in the asymptotics of certain expectation values of the Heisenberg- 
picture operators e'^Nt/h^{p)^-iHt,t/h^ as N ^ oo. If Hn is chosen as in (1.1), (1.2), and 
A^^^ is chosen as in (1.4), the limit — > oo is the usual mean-field limit; see [He, Sp]. 

Our first main result is the following 



Theorem 1.1 . Let h> and t > be fixed, and let w G L°^(R'^). If'i>j^fi{xi, . . . ,xj^) = 
^l^{xl) ■ ■ ■ '(^{xn) is a normalized "coherent" (i.e., product) initial state, then 

lim {^n,Oj 

lim {^N,t,AP^N,t) = (^p,t,a(P)M/p,t) =: a^P\i;t) 

Here'^N.t is again a coherent state, i.e., N,t{xi, . . . , xn) = '4't{xi) • • • iptixN), and '^p,t = 
^ N=p,t, where ipt is a solution of the Hartree equation 

ihdt^t = -^^i^t + {w*\i^t?)^t (1.6) 

with initial condition ipt=Q = V'- 
Remarks 

1. For large N , the quantum evolution e~'^^^^^^'^ can be replaced by the nonlinear 
single-particle evolution ^'Ar^t(xi, . . . ,xn)- Particle interaction effects are translated 
into the nonlinearity of this evolution. This justifies interpreting the limit ^ oo as 
a mean-field limit. 

2. The corrections to the limit in (1.5) are 0{1/N). 

3. Since lim N{N — I) • • • {N — p + 1)/NP = 1, the second equality in (1.5) follows 

N—>OD 

easily from (1.4), because 



4. Theorem 1.1 was first proven in [He], see also [GiVe]. A new proof was given in [Sp] 
and extended to more general classes of two-body potentials, including the Coulomb 
potential, in [EY], [BGM], [BEGMY]. The proof in our paper is quite different. It is 
inspired by a second-quantization formalism to be published elsewhere. It enables us 
to tackle the problem of obtaining convergence estimates uniform in Planck's constant 
h, as we now proceed to discuss. 
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It is well known that, for Wjy as in (1.2), the classical dynamics of N particles tends to 
the dynamics defined by the Vlasov equation, in the limit N ^ oo. More precisely, if pat 
denotes the empirical distribution, namely 

1 ^ 

1=1 

where . . . ,xj\[{t); ^i(t), . . . , (,N{t)) is a solution of the classical equations of motion, 

then, in the limit N — > oo, pj\f tends weakly to ft{x,(,)dxd(^, where ft{x,^) is a solution of 
the Vlasov equation: 

dtft = -^-VJt + VxVeffV^ft (1.7) 
Veff{x,t) = I w{x-y)ft{y,C)dydC , (1.8) 

see [BH] . It is natural to ask whether this convergence result is related to that of Theorem 
1.1. Our next result provides, under very restrictive assumptions on the two-body interac- 
tions, a partial answer to this question. First, we define a restricted class of interactions. 
For 0" > 0, we define the spaces 

Ll^p := {/GLi(M6p)|e-N/eLi(M6P)}, (1.9) 
Aa,p := {f eLHR^P)\e''^'^f(^LHR^P)}, (1.10) 

Here Xj £ M^, G M^, j = 1, . . . ,p, and 

z := iXp,Ep) G ]R3p ^ ]^3p. -.= (xi, . . .,Xp), Ep := (6, . . .,Cp), 

p 

:=I^(kjl + l6-l); 

/(s), s := {S, S) G M^p X M.3p ^.j^g Fourier transform of /. 

We further denote by <I>^ : {Xj\[;En) •— > (XAr(t); HAr(t)) the flow generated by Hf^, 
where Hf^ is the classical Hamilton function corresponding to the operator H]\f. 

Definition 1.1 We define by: 

1. 

Wl^{XM,EN;t)= (1.11) 

(27r)-3^ / e^<^^-^>^iv(XAr + hYN/2,t)^N{XN - hYN/2,t)dYN 
the Wigner distribution of the N -particle normalized wave function ^ n{Xn 
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2. 

= / T^*^(X^,S^;t)dX,v-i'iHjv-i, (1.12) 

yR3(]v-i) 

the j— particle Wigner function ({N — j) -marginal distribution of the N -particle Wig- 
ner distribution). 

3. 

Wmx, t) = (2^)-3 f Sy^^Htix + hy/2Wt{x - hy/2) dy, (1.13) 
the Wigner distribution of the solution ipt{x) of the Hartree equation. 
Our second main result is 

Theorem 1.2 . Let w E for some a > 0. Let ^n,o be a product state. Set e := 

||i(;||ooi. Then, for fixed p, there is a constant Cp > independent of h such that, as an 
equality between tempered distributions, 

W^-{X„E,;t) = fl Wmxj,^,;t) + § + O (e^Vv^) , (1.14) 

as N ^ oo. 
Remarks 

1. It is known that W{ip){x, ^; t) converges in 5'(M^) to a solution ft{x, of the Vlasov 
equation, as h ^ [NS]. It is also known that 

p 

hm I<^(Xp, Hp; i) = n fti^j^^j) 
whenever N ^ oo entails — > 0, as in the case of the Kac potentials [NS],[GMP]. 

2. Result (1.14) shows that, up to an exponentially small error independent of h, the 
mean-field convergence towards a single-particle nonlinear dynamics holds uniformly 
in h. 

3. The classical limit is equivalent to the limit of heavy particles. We set h = 1 in (1.2), 
but let the particle mass m become large. We impose the condition that the kinetic 
energy per particle be independent of m, namely mvf = 0(1), i.e., \vi\ = 0{l/ y/rn), 
for all i. This suggests to rescale time as t = ^JraT. Then the Schrodinger equation 
becomes 

. ^ A 1 ^ 

^ i=i ij=l 

which is equivalent to (1.1)-(1.3), for h = 1/y/m. 
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2 The ^ oo limit: convergence estimates 

2.1 Kinematical algebra of " observables" 

The above systems can be described by a kinematical algebra of operators, the quantum 
mechanical analogue of the algebra of functions on phase space of a classical system. 

Let n^P^ := LliR^P), < p < N, N eN. Let a^P^ be a bounded operator on n^P\ and 
a^P\xi, . . . , Xp; yi, . . . , Up) := a^P\Xp\ Yp) be the tempered distribution kernel in S'i^P x 
M?P) associated to a^^^ by the nuclear theorem: 

{a^P^^^P^){Xp)= [ a^P\Xp-Yp)^^P\Yp)dYp (2.1) 

where f'^\Yp) G LI{W^p). Then [a^P^)* has the distribution kernel a^P){Yp;Xp). 
To we associate the operator on L|(]R'^^) specified in (1.4). Explicitly: 

4)(a(-))^(^) := f^) / K{{.lY^)^^^){Y^)dY^, (2.2) 

with 

K(x^,y^) = a(p)(Xp;yp)(^(x^_p - yjv-p), ^^^\Yn) g l|(]r3^) (2.3) 

If o(P) is bounded on ^(p) then A'^\a^P'>) is bounded on Ti^^l Since IIP5II = 1 and 

— — < 1, we have that 
\p J NP - 

< / ( [ \v^^''HYp;X^.p)\^dYp) dX^.p 

We set 

i:jv := (AS^)(a(f))|a(P) e S(w(p)), p = 0, 1, 2, . . .) C B{n^^^). (2.4) 
The following statement is easily verified. 

Proposition 2.1 . The map a^P'^ ^ A^^\a'^P'^) is linear, and {A'§\a^P'^))* = A^^\{a^P'>)*), 

2.2 The Schwinger-Dyson expansion 

Given a bounded operator A^^^ acting on 7{^^\ p < N,we denote by A^l^ the corresponding 
Heisenberg-picture operator with respect to the free time evolution e^^'^^^^, i.e., 

^ ^iHo^t/n^{p)^~iH"^t/n (2.5) 
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We further denote by A^^\t) the corresponding operator with H'^ replaced by H]\f, namely 

:= eiHMt/h^{p)^-iH^t/H^ (2.6) 
and by A^f^^{t,s) the two-parameter operator family 

s) := e'H^'/''^-iH%t/n^(j^)^^iH%t/n^-iH^t/n (2.7) 

Then we obviously have 

A^\t) = A%{t,s) (2.8) 

' s=t 

We denote Wjy by W. Iterating the identity: 

Af^^it, s) = A% + 1 f e'"^'-l^e-'<'^'^[Wt, , ^f^Je^^V^g-i^^iv//^ dt, (2.9) 
^ Jo 

we get that 

C50 / ' \ ^ ft ft _1 

Af},{t,s) = A% + Yl(^^) I dh...J^ ' dt,,[Wt„,...,[Wt„A%]...] (2.10) 
and finally, setting s = t, we obtain the Schwinger-Dyson expansion 



n=l 



A^N\t) = A% + f2[l)) / / ■■■ r'\wt^,...,[Wt„Al%]...]dtn...dh (2.11) 



From now on, we drop the index N in the Heisenberg-picture operators with respect to the 

(p) (p) 

free evolution, i.e. we use the abbreviation: A^ '■= Al . 

(p) 

The boundedness of Af^ and of the interactions Wt^ implies the boundedness of all multiple 
commutators, with 

• • • , m^A^f^] . . .]||^(^.) < {2\\w\\^,^,/hr\\A%\\^,^, 

<(2||Ty||„w//i)"||a(P)||^(p), 

for A%^ = A%\a^P'^). By (1.2), ||W^||^(iv) oc A^. Hence, for fixed N and /i, the series is 

norm-convergent, for all t > 0. The time integrations yield a factor — , so that the norm 

n! 

of the series in (2.11) is bounded by exp [2||VF||-^{jv) • ||-^(p) • 

These estimates are obviously not adequate to investigate the ^ oo or the ?i — > 
limit, let alone to prove uniformity in h. 
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2.3 The N ^ oo limit 

We exploit the structure of the commutators on the right-hand side of (2.11), the symmetry 
of wave functions in LKM'^^), and the fact that each term in only acts on p arguments 
of a wave function, so that many commutators will vanish. Note that 

TV JV 
i<j i<j 

where 



Therefore 



1 ^ 



i<j 

p N 1 P 



1=1 j=p+l i<j 



^E^!^'''^(K^+^«^"^]) + ^E4^(K,«^^^]) (2-14) 



In more precise terms, the expression 

[Ws,A^^^] = ^E4^'^([^r\a(^)]) + ^E4HK,«^^^]), (2-15) 

i=l i<j 

holds as an operator identity on TL^^\ In second-quantization language, the first sum on 
(2.15) corresponds to tree graphs, the second one to loop graphs. Next, we insert (2.15) in 
(2.9) and perform a second step, but only for the first sum in (2.15), leaving the second 
one unchanged. To keep our notation compact, it is useful to introduce the notion of tree 
amplitudes of n— th order, recursively defined in the following way 

,(0.)=,(P). ,a„ = l'£\u;;r,,fc^ej, n>l (2.16) 
Then expression (2.15) becomes 

The first term is 0(1), while the second one is of order p{p — 1)/N (for fixed h) and 
is therefore suppressed by a factor 1/N. Performing (k — 1) iterations only for the tree 
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amplitudes, we conclude that 

k 

N 



k 



n=l 



*""4^"^(5aj^in,-..rftl, (2.18) 

>i=i -^0 ^ 

(2.19) 

^iH%t„/h^~iHr,U/h-^l^j^^ p. j^yHNt„/h^-iH%t,,/h _ ^ ^ ^ ^2.20) 



i/(.,p;iV):= Y: 4^^"'^(K,ff£L_J) (2-21) 

Equation (2.17) is most easily verified as follows (think of A^"^" ^\9^i~^'^„-i') a, p+n — 1- 
particle operator replacing A^'lf in (2.15) ): 

. [Wt^ " (C..1-1 )] {N - {p + n - 1)) ."Is^ 

/i(p+")/'r *p+" iiP 



where 












p+n-1 i 




^t,N 


rt f-tk 
Jo Jo 


with 





h N 

i=l 

p+n-1 



J>i=l 

p+n— 1 



j>i=i 

2.4 Control of the expansion, small time, h fixed 

First, we prove a bound on the norm of A^^~^"'\gl'^'^\^J 

PS^^"^(^?ajll>.(-) < liP + n~l)\M^ ll4^"-'^(5fc.!?„JII„(iv, (2.22) 

This follows from the unitarity of the free time evolution and the boundedness of the 
interactions, = HujUoo- The bound (2.22) then yields recursively 

PS?^"^(ffKjlk(-) < (f + ^ - l)(f + n - 2) • • • (p + l)p Qlklloo)" Il4^||^(iv) 

- ^^7r^Qll^ll-)"ll«^"^ll«(.)> (2-23) 
independently of all time indices. 



Considering the expansion (2.17), we have that 

- 

lo Jo 



n=l 
oo 



< 



n=l ^' ^' 



n=l 



W no t 



(2.24) 



because ; — ; — < 2"""*"^. The series on the R.S. of (2.24) converges for \t\ < [ —\\w\\, 

p'.ni \h 

The third term in (2.17) is bounded similarly. Let 



Then 



oo p+n—l j- -f- _j 

n=l i<j=l -"^ •'^ 



< 



N 



Jo 



1 ^ {p + n-l)^ 2 
2 W 

n=l 
^ oo 

^]7E(?'+'^-i)'(^'+^-2)- 

n=l 



•P 1 ^ll'f^lloo 



n-1 



^l|ooP(^)|l7,(iV)^< 



' n=l 



N 



2^E(fiii"^ii-i*ir 



Therefore 



n=l 



- ■ n=l 

The remainder term in (2.17) clearly vanishes, as A; ^ oo. To summarize, we have proven 
the following result. 



Proposition 2.2 . Let \t\ < 



^ * n=Jo Jo 



2.5 Convergence for all times, h fixed 



4^"^(5&ijdtn...dtl+0(l/iV) 



(2.26) 



We assume that the statement of Theorem 1.1 holds up to some time T independent of p, 
i.e., 

hm e^^'-'^/'^^S^)e-^^^^/^M/;v,o)H(^) = a^'^HM (2.27) 



7V-»oo 



Let us proceed one step further in time, with t < (4||i(;||oo/?i) ^- On account of (2.17), we 
have that 

^iHN{T+t)/hj^(.p) ^-iHN{T+t)/h ^ 

^iHMT/h^-iHNt/hj^{p)^-iHNt/h^-iHNT/h ^ ^iHMT/h^(p)^^-^^-iHMT/h (2.28) 

^ ft 

+ E / • • • / e^^-^/^Ai^+")(5K,Je-^-^/^dt. ...dt, + 0{1/N) 

n=l 



This expansion is norm convergent, by (2.17) and the unitarity of g^^NT/h^ Taking expec- 
tation values in ^'at^o we get, as above, 

1 Jo Jo 



n=l 

+ 0{1/N) 



Hence, by the inductive assumption and the norm convergence of the series 

Af^oo 

X: lim / . . . / (^^,, e^^-^/'^^g+")(<7£St Je-^'''^^/'^^,*)7.(-) dtn... dh 



n=l 



We postpone to Section 3, below, the proof that actually 



lim ... {^N,uA%'^''\gT,ZJ'^M^t)n^N)dtn...dh (2.29) 

^^^t^iJo Jo 



•t+T) 



and this ensures that the convergence is global in time. 

2.6 Control of the expansion, uniformity with respect to h 

Given a symbol t{x, ^) E Aa,p, we denote by T the corresponding Weyl operator. Its action 
on vectors ip £ 5(M'^^) is given by 

(rV')(x) = -^ / / r[(x + y)/2,e]e^«"-^)'«)/^V(y)« (2-30) 

In general, T is a semiclassical pseudodifferential operator. Let us recall some relevant 
results (see e.g.[Ro]). 
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1. If r G L^(M'^P X R^P) then T extends to a continuous operator on L'^{M.^p) with 
\\T\\ < \\t\\x^i; hence ||T||2,2^j;^2 < \\T\\^ := ||t||o-, where r is the Fourier transform of 
r, and 

||r|U := [ \f{s)\e''^'^ ds. 

Obviously, ||r||/,i < ||t||o-. 

2. U w £ Aa,p,g G Aa,p, then the symbol of the commutator [W,G]/ih is the Moyal 
bracket {w,g}M- Hence the multiple commutator [W, [W, . . . [W,G] . . has as 
its symbol the multiple Moyal bracket {w, {w, . . . ,{w,g}M ■ ■ ■}m- 

We recall that, given g,g' G Aa,p, their Moyal bracket {g,g'}M is defined as 

{g,g'}M = g#g' - g'#g, 

where ^ is the composition of Weyl symbols. In the Fourier transform representation 
the explicit expression of the Moyal bracket is (see e.g. [Fo],§3.4): 

{{g,g'}M)''{s) = l [ g{s')g' {s-s') sin [h{s-s') A s'/2]ds\ (2.31) 

where, given two vectors s = (vjw) and = {v^,w^), s A := {'w,vi) — {v,wi). 

3. If the observable T has symbol t{x, (^), then the Heisenberg observable Tf has symbol 
(r o <I)^)(x, ^). Here ^f{x,^) = (x + .^t,^) is the free flow with initial conditions (x,^). 
In particular, (r o $^)(x,^) G Acr,p whenever r € Aa,p- 

Under the present assumptions, it can be proven, starting from the expression (2.31) (see 
[BGP], Lemma 3.2 ), that the following estimate on the Moyal bracket holds: 

\\{w,g}M\\a^5 < ^TI^WMl^rhWa, 0<6<a (2.32) 

For the convenience of the reader we reproduce here the proof of (2.32). Since {s — s^) As^ = 
s A s^, and |s A S"*^! < |s| • \s^, by definition of the ^o--norm and (2.31) we get: 

\\i{g,g'}Ai)\\a-5 = 

If e("-^)l^lds / \g{s^)g' (s-s^) sin {h{s-s^) A s^)/2\ds^ 

f / e("-^^(l^l+l"'l)|5(s)5'(si)sin(n(s Asi)/2|(isi 
</ e(^-^)l^l|5(s)|(is / e^''~'^\'\'{s')sAs'\ds^ 

Jr'^p Jr'^p 

< [ e("-^)l''l|5(s)||s|ds / e^''-^^^''^\g'{s^)s\\s^\\ds^ 

jRSp Jr6p 
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whence the assertion because xe ''^ < — , Vx > 0, V(5 > 0. 

ed 

Let now Qj. := {gr-i,w}M,r > l;gi = {g,w}M- Then, applying (2.32) r times, we can 
write: 

WdrWa^rS < (^-^^ \\w\\l Mia (2-33) 

These results immediately yield the following bound. 

Lemma 2.1 . Let the operator a^^^ be the Weyl quantization of a symbol Ta{x,^) G Acr,p 
for some a > 0. Then there is L{p) > independent of h such that 

P!?^"^(5SJII^(-) < L-n}!^{2M\^T\\a^^\ (2.34) 

Proof 

Denote by Q'flf\.tn symbol of g^''^}tn- Using definition (2.16) and the estimate (2.32) 
we get the uniform estimate corresponding to (2.23): 

Wri'^^P) II ^ 2(p + n— 1).. .. ||^(n,-l,p) 11 n ^ r ^ riK\ 

WtMLtr.-,Nh~^5n < ^272 ll^ll<^ ll^i,il,...,tn-lll'^' < (5„ < (J (2.35) 

The recursive definition (2.16) allows us to use the recursive estimate (2.33). We get 

ii<";l,tJi-n^» < neXr^^^^MUa^^^L (2.36) 

Setting 6n '■= — we get the bound (2.34) on account of the majorizations 
2n 

II ^(p+n) , (n;p) nm ^ ||^(n,p) n n n II II ^ II II 

PtV (9ti,...,t„)ll>^(iV) < ||^t,ti,...,t„IU/2, \\w\\l2^l^ = \\w\\oo < \\w\\„. 

This proves the Lemma. 
Remark 

The uniform control in h introduces an extra n!^ divergence with respect to the fixed-^ 
estimate (2.23). 

We now obtain uniform estimates of the three terms in expansion (2.17). 

Lemma 2.2 . There exist constants Mi > 0, M2 > 0, > 0, Li > 0, L2 > 0, L3 > 0, 
independent of (h, t) and N, such that 

k 

< Mi||a(P)|Uj;L^n!2(||^||^t)" (2.37) 

n=l 

\\Qt)t\\nW < MsWa'^P^lU^nlWM^tr (2.38) 
ll^i'llwciv) < MsWa^P^W^LlklWM^t)'^ (2.39) 
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Proof 

Inserting the estimate (2.34) in the expressions (2.18,2.19, 2.20) we get, on account of 
unitarity of Uo(t): 



k ft ft 



n-l 



< ||a(P)|UV(2L||«;|U)"n!3 /.../■"' 

k 

< ||a(f)|U j;(2L||^«|U|t|)"n!2. 

n=l 

The last inequahty comes from performing the time integrations, which are majorized by 
a factor in (2.37,2.38) and by a factor Itj'^/A;! in (2.39) (proven with the help of the 

same argument). This proves the lemma. 

Using this result, we can easily prove the uniform version of the expansion (2.17). 

Proposition 2.3 . Let e := || t(7||ooi' Then, in the scime assumption of Lemma 2.1 on the 
operator there exists k = k{e), A = A(e) such that 

^^H^t/n^{p)^-iH^tn ^ ^(rt +5(p).fc + ^(p);'= + A, (2.40) 



where 

bI^n' = E / • • • / ^N^""^ )dtr.--- ■ (2.41) 



n=l 



Jo 



Here B^^^ fulfills the majorization (2.37), and 

||i2lj''|l7^(iv) <M3e-^^/v^ (2.42) 

Proof 

The estimate (2.38) and a standard Nekhoroshev-type argument show that the choice 

k{e) := ^ = (2.43) 

minimizes the divergence of R[^1}^ ■ A straightforward computation then yields (2.42). By 
definition of Q^}^^ we get the uniform version of the estimate (2.21), whence 

k{e) 

A(e) :=p2P\\a'-P^\\^^n]\2e)''' < p2P€^^/\e/2y^/^' 

n=l 
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3 Connection with the Hartree equation and proof of the 
theorems 

We wish to prove that the representation of the evolution obtained in Proposition 2.2 
coincides with the evolution generated by the Hartree equation in the limit N ^ oo . 
For this purpose, we recall that the Hartree equation is Hamiltonian. We define the func- 
tional 

^(V',V^) = -v/ |VV(^)pdx + W(^,V^), (3.1) 
2 Jrs 

for ip G H^{R^), where 

y^{'^,'^) = \( Tp{x)'il:{y)w{x - y)'il:{x)ij){y) dxdy (3.2) 

If ip{x),'ip{y) are considered as canonical variables with Poisson brackets 

{tp{x),tp{y))} = ih6{x - y), {ip{x),ip{y)} = {V'(x), V'(y)} = 0, 

then (3.1) is the Hamiltonian functional generating a time evolution of functionals on phase 
space equivalent to the Hartree equation. Namely, if A{ijj) is a functional and At denotes 
its time evolution, one has that 

dtAtW = ^{n,At}W 

Choosing A = {(t>,ip), 4> ^ C'o"(I^'^)) then At{il^) = A{ilJt)i where il^t is a solution of the 
Hartree equation 

ihdt^t = -^Mt + {w*\M^)^t (3.3) 

Define the free flow ^>?(^) := At of A by 

At = ^(e^^*/V) 

and denote by ^t{A) the interacting flow. Formally, the interacting flow is given by the Lie 
expansion in the interaction representation (analogous to the Schwinger-Dyson expansion 
of Section 2.2). Indeed we have the following result: 



Lemma 3.3 . ^t{-A) admits the formal expansion 

^t{A) = A + ff; f^... f\m^ ...{Wt,,A}.. .} dtn... dt^ (3.4) 
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Proof. To see this, we consider the dynamics in the interaction picture. We set 

At := o $° t(^) 

Then 

n 

where Wt := <I>j (W) is the free evolution of W. After integrating in time we get 







whence 



Iterating this identity, we obtain the series (3.4), and this concludes the proof of the Lemma. 
The desired identification is based on the following proposition 

Proposition 3.4 . Let G H^(R^), and let be a product state, i.e. 

I 



Xs) 
s=l 



Then, for all N > p, 



(3.5) 



^) {Wt„...{Wt,,A}...m, 



where 

r P 

A = a^P\il;) := / ij{xi) ■ ■ ■ ip{xp)a'^\xi, . . . ,Xp]yi, . . . ,yp)i){yi) ■ ■■tp{yp) dx^dy^ 

k=i 

Proof. We have that 

o ^\{{WuA]/h)l^^ = {W,A}/h = dtA\t=o 

Define the projection p := \iIj){tP\. Then, denoting pt := ° ^-t{p), we have that, for all 
n > 1, 

At = TviApfn-, dtA\t=o = dtTr{Apf^)\t=o. 
Therefore, in the interaction picture 

dtpt = TripfW^' - Wrpf)/h 
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and in the same way we get 

1 " 

1=1 

It is then easy to check that 



and this concludes the proof of the Proposition. 

We are now in a position to prove our main results. 
Proof of Theorem 1.1 

Consider the expectation value of the expansion (2.26) in a coherent (i.e., product) state: 

••• / (*^,AS+"^(9Kj^;vW)dt....dti + 0(l/iV). 
JO Jo 



n=l 

By definition of ^'tv, 



N{N -1)...{N -p+l) 



n=l 



Jo 



+ Y, ■■■ i^n+p, 9?, ,tn^«+p)w("+rt dtn ■■■dti+ 0{1/N). 



Since the series is norm- convergent, the limits N ^ oo and n ^ oo can be interchanged. 
Then 



hm {^N, e^^^*/^A^e-^^^*/^^jv)w{^) = (^p, af^^p)^,,, + (3.6) 

I — >O0 

1 Jo Jo 



n=l 



where the last equality follows from formula (3.5). 
Proof of Theorem 1.2 

If, instead of (2.26), the representation (2.40) is considered, the above argument yields 

(3.7) 



16 



Given any bounded operator A on L^(IR^') with (Weyl) symbol a^ix,^) : S{ 
where S is the Schwartz space of rapidly decreasing functions, its matrix elements can be 
expressed in terms of the symbol and of the Wigner function by the following well known 
formula (see e.g.[Fo]): 

{^,A^)=[ aA{x,OW^{x,Odxd^ (3-8) 
where W<i/{x,S^) is the Wigner function of the state Therefore, in our case 

iV(iV - ij . . . (7V -p+ ij J]g3pxR3p 

where W^"" {Xp,Ep,t) is the Wi gner function corresponding to the time evolution, 
Q^^Nt/h^ ^ ^ of the product state ^'at^o = i^ixi) ■ ■ ■ tp{xj\f). The N—p variables {X]\f-p, Hat^p) 
are integrated out. By (3.7) and (1.13), we can take the N ^ oo limit and write 

/ aA{Xp, Ep)W^^{Xp, Ep, t) dXpdEp = 

/ (JA{Xp, Ep) TT W^{xi,Ci; t) dXpdEp + 0(e-i/v^). 

Since this formula holds for any aA{Xp, Ep) G 5(M'^p x R'^^) n Aa,p, the assertion is proved. 
Proof of formula (2.29). 
By (3.5), we have that 



N) 



{Wi„...{Wt,,^}...}(VT) 

which yields formula (2.29), by Lemma 3.3, on account of the uniform convergence of the 
series. 
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